On the elastic energy and stress correlation in the contact between elastic solids with 

randomly rough surfaces 
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When two elastic solids with randomly rough surfaces are brought in contact, a very inhomoge- 
neous stress distribution (t(x) will occur at the interface. Here I study the elastic energy and the 
correlation function (cr(q)(T(— q)), where (T(q) is the Fourier transform of o"(x) and where (..) stands 
for ensemble average. I relate (a"(q)(T(— q)) to the elastic energy stored at the interface, and I show 



that for self affine fractal surfaces, quite generally ((T(q)(j(— q)) 
exponent of the self-affine fractal surface. 



-(l+H) 



where H is the Hurst 



The contact between elastic solids with surfaces with 
roughness on many length scales is a topic of great prac- 
tical importance. Particularly important are the area of 
real contact A, which in most cases determines the slid- 
ing friction force, and the asperity induced elastic energy 
Uci stored at the interface, which is directly related to 
the interfacial separation [1]. 

In a recent paper Campana, Miiser and RobbinsQ 
have presented numerical simulation results for the stress 
distribution at the interface between two elastic solids 
with randomly rough (self affine fractal) surfaces. They 
have shown that the stress correlation function 

(a(q)a(-q)) ^ q-(^'+"^ 

where H is the Hurst exponent of the self afhne fractal 
surface. Here I show that this is also the prediction of 
the contact mechanics theory of Persson[3j. I also derive 
some other results related to the elastic energy and the 
accuracy of the Persson theory. 

Consider the contact between two elastic solids with 
rough surfaces, but which appear flat at low resolution. 
We can write the elastic energy stored in the vicinity of 
the asperity contact regions as 



Uei = \ I d'x (a(x)u(x)) 



(1) 



where m(x) and ct(x) are the normal displacement and 
the normal stress, respectively. We write 



a(x) = J (fq a(q)e'i-^ 



(2) 



and similar for u(x). Substituting this in (1) gives 



[/ei = (2^)2- / d^q (a(qV(-q)) (3) 



and 



C/el = (27r) 



For complete contact u(q) = 
tion 

(/l(q)/i(-q)) 
gives for complete contact 
t^cl — : — 



d?q q{u{ci)u{-Ci)) (6) 
= /i(q) and using the defini- 

^0 



Ml) 



d?q qC{q) 



Now, for incomplete contact I define W{q) so that 



d\ qC{q)W{q) 



(7) 



(8) 



(9) 



I have argued elsewhere^ that W{q) = P{q) is the rela- 
tive contact area when the interface is studied at the mag- 
nification C — q/qo- The qualitative explanation is that 
the solids will deform mainly in the regions where they 
make contact and most of the elastic energy will arise 
from the contact regions. In particular, for complete con- 
tact P{q) = 1 and in this limit using W = P = 1 would 
be exact. I will show below that W{q) ~ P{q) follows 
directly from the contact mechanics theory of Persson, 
but for the moment we simply define W{q) via (9). 
Comparing (5) and (9) gives 

(a(q)a(-q)) = (^^J Aoq'C{q)Wiq) (10) 

In what follows we need 

d'q (a(q)a(-q)) 



Next, using that[ 



a(q) = -E*quici) 



gives 



(4) 



f/e, = (2^)2-i- I d\ g-i(a(q)a(-q)) (5) 
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(27r)4 



(27r) 



'\\p-i^-i:^-^') 



d'xd^x' / d'q (a(x)a(x'))e 



d^x (a2(x)) = -^(a^) (11) 



(27r)2 



Substituting (10) in (11) gives 



where 



« 9 (rr'^) 



(12) 
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All the equations (1)-(12) presented above are exact. 

Now, it is easy to calculate (cr^) approximately using 
the Persson contact mechanics theory . The basic equa- 
tion for the stress distribution P(cr, C) 



where 



/(C) = jE*'q,q^Ciq) 



(13) 



(14) 



where q = goC- Multiply (13) with and integrate over 
a: 



/(C) 



= 2/(0/ daP{a,0 
Jo 



(15) 



where I have performed two partial integrations. Now, 
note that 



Jo ^0 



(16) 



which I also denote as P{q) {q = qoC) for simplicity. Us- 
ing the definition of f{() we get 



(^') = ^E*' / dq q'^C{q)P{q) (17) 



Substituting (17) in (12) gives 



dq q'Ciq)Wiq) = / dq q'C{q)P{q) (18) 



Thus we get W{q) = P{q)- The resulting equation for the 
elastic energy was used in Ref. h\, but only qualitative 
arguments was given for its validity (see above). Here we 
have proved that it follows rigorously within the Persson 
contact mechanics theory. Substituting this result in (10) 
gives 

(a(q)a(-q)) = [^^^ ' A,q^C{q)P{q) (19) 

Assume now a self affine fractal surface. In this case 

C{q) ^ g-2(«+l) (20) 
If we assume that the relative contact area is small then 
P{q) « [7rG(q)]-V2 (21) 



dq q'C{q) 
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2{1-H) 



2{1-H) 



when q >> go- Thus, 



P{q) ^ q 



H-1 



(21) 



(22) 



when the contact is small and q » go- Substituting 
these results into (19) gives the q-dependence 



(a(q)a(-q)) 



(23) 



in good agreement with recent numerical studies of 
contact between elastic solids with randomly rough 
surfaces 0|- 

Detailed analysis of Molecular Dynamics Finite 
Element Method [6j and Green's function Molecular 
Dynamics Q calculations indicate some small deviations 
from the prediction of the theory of Persson. Thus, the 
contact area calculated as a function of the squeezing 
pressure tend to be about ^ 20% larger at small pressures 
(where A varies linearly with p) than predicted by the 
Persson theory. When the squeezing pressure increases, 
the deviation becomes smaller and vanish at complete 
contact. Similarly, the elastic energy in the asperity con- 
tact regions may be overestimated in the theory of Pers- 
son at low squeezing pressure (see i), while the differ- 
ence decreases at higher squeezing pressures and vanishes 
for complete contact. Here we will use the theory above 
to show that these two facts are in fact related. 

We focus on small contact where the area A of contact 
depends linearly on the nominal squeezing pressure p. 
Since Persson's theory predict smaller contact area than 
numerical simulations, the pressure in the contact regions 
will be higher and we expect that the prediction of the 
Persson theory for (ct^) will be higher than observed in 
the numerical simulations. In fact one can show that 



{a') « {a')pAp/A 



(24) 



where the index "P" refer to the Persson theory. Thus, 
using (12), which holds exactly, we conclude that W{q) w 
P{q)Ap/A. Thus, at low squeezing pressures, using the 
Persson prediction W{q) = P{q) will overestimate the 
elastic energy by a factor of Ap/A, i.e., by about ~ 20%. 
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